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A theoretical model for driving a two qubit system to a stable long-lived entanglement is discussed. 
The entire system is represented by two atoms, initially in ground states and disentangled, each one 
coupled to a separate cavity with the cavities connected by a fiber. The cavities and fiber exchange 
energy with their individual thermal environments. Under these conditions, we apply the theory of 
microscopic master equation developed for the dynamics of the open quantum system. Deriving the 
density operator of the two-qubit system we found that stable long-lived quantum correlations are 
generated in the presence of thermal excitation of the environments. To the best of our knowledge, 
there is no a similar effect observed in a quantum open system described by a generalized microscopic 
master equation in the approximation of the cavity quantum electrodynamics. 
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I. INTRODUCTION 

Entanglement (verschriinkung) introduced in physics 
originally by Schrodinger [l| and considered a native fea- 
ture of the quantum world, is the most outstanding and 
studied phenomenon to test the fundamentals of quan- 
tum mechanics, as well as an essential engineering tool 
for the quantum communications. However, entangle- 
ment is a property that is hard to reach technologically 
and even when achieved, it is a very unstable quantum 
state, vulnerable under the effects of decoherence, any 
dissipative process as a result of the coupling to environ- 
ment. Conventionally these effects are considered mainly 
destructive for entanglement, nevertheless some recent 
studies of this subject attest results different from the 
common conviction, even appearing as counterintuitive 
at first glance 



An alternative approach to measure the entire correla- 
tions in a quantum system was suggested in Refs. 0, Q. 
For example, by using the concepts of mutual informa- 
tion and quantum discord (QD) the quantum correlations 
may be distinguished from the classical ones. Further the 
QD could be compared to the entanglement of formation 
(E) in order to find if the system is in a quantum in- 
separable state (entangled), or in a separable state with 
quantum correlations, sucha as QD |8l-lll|. Such an anal- 
ysis is considered in this paper. 

The inclusion of the interaction of the system with the 
environment plays an important role in physics, implying 
a more realistic picture because the dissipation is always 
present in the real devices. In the proposed study we 
deal with atoms, cavities and a fiber in the framework 
of the physical model suggested in Ref. [T^], which at- 
tracted a high interest for quantum information appli- 
cations and subsequently discussed details from different 
aspects [l3l - [l6j . As a basic model, we consider the one 
recently analyzed in Ref 



tions for a very special case, i.e 



17| and extend the calcula- 
when the atoms are 



initially disentangled and in the ground states while the 
fields are in vacuum states and coupled to the reservoirs 
at finite temperatures. The entire system is considered 
open because of the leakage of the electromagnetic field 
from the cavities and fiber into their own thermal baths. 
Therefore, we ask ourselves the following question: Is it 
possible to generate atomic quantum correlations by the 
processes of absorption and exchanging excitations with 
the thermal reservoirs? In the following we present the 
model and detailed analysis in search of an answer. 



II. MODEL 

We present here the model schematically shown in 
Fig m and recall the basic equations that lead us to the 
effect we are looking for. Hence, one considers two qubits 
(two-level atoms) interacting with two different and dis- 
tant cavities, coupled by a transmission line (e.g., fiber, 
waveguide). For simplicity we consider the short fiber 
limit: only one mode of the fiber interacts with the cav- 
ity modes [T5| . 

Now, let us define a given state of the whole system 
by using the notation: \i) = \Ai) ® |^2) ® \Ci) ® IC2) ® 
\F) = \A1A2C1C2F) , where ^^=1,2 correspond to the 
atomic states, that can be e{g) for excited(ground) state, 
while Cj=i^2 are the cavity states, and F corresponds to 
the state of the fiber. Both Cj=i^2 and F describe a 
or 1 photon state. The Hamiltonian of the composite 
system under the rotating-wave approximation (RWA) 
reads (with h—1) 
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+ {djS+aj + i^asa] + H.c}j , (1) 
where 03 is the boson operator defining the fiber mode, 
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FIG. 1: Two atoms trapped in distant coupled cavities. The 
cavities and transmission line exchange the energy at the rates 
7i, 72 and 73 with their baths having the temperatures Ti, 
T2 and Ti, respectively. 



01(02) is the boson operator for the cavity 1(2); wq and 
uja are the cavity and the atomic (fiber as weW) frequen- 
cies, respectively; gj [v) the atom-cavity (fiber-cavity) 
coupling constants; and Sz, are the usual atomic in- 
version and ladder operators, respectively. The model 
is studied under the assumption of a single excitation 
in the system of atoms and fields, and using the above 
mentioned notation, the state-basis of the system be- 
comes: |1) = legOOO) , |2) = \gglOO) , |3) = |.g,g001) , |4) = 
|.g.g010) , |5) = \geOOO) , |6) = \ggOOO), where the last vec- 
tor is required by the existence of the excitation's leakage 
to the reservoirs. Hence, it is straightforward to bring the 
Hamiltonian Hg in Eq. ([T|) to a matrix representation in 
the state-basis [l7j . 

To simulate the dynamics of the given system, one con- 
siders the approach of the microscopic master equation 
(MME), developed in Refs. [H, [H in order to describe 
the system-reservoir interactions by a Markovian mas- 
ter equation. This description considers jumps between 
eigenstates of the system Hamiltonian rather than the 
eigenstates of the field- free subsystems, which is the case 
in many approaches employed in quantum optics. There- 
fore, we assume that the system of interest, i.e. the 
atoms, cavities and fiber are parts of a larger system, 
composed by a collection of quantum harmonic oscilla- 
tors in thermal equilibrium. The external environment 
represents the part of the entire closed system other than 
the system of interest. Between each element of the sys- 
tem and its own bath one may identify different kind of 
dissipation channels. In cavity quantum electrodynamics 
(CQED) the main source of dissipation originates from 
the leakage of the cavity photons due to the imperfect 
reflectivity of the cavity mirrors. A second source of dis- 
sipation corresponds to the spontaneous emission of pho- 
tons by the atom, however this kind of loss we consider 
small and neglect in the model. Following the common 
procedures [la, [l^, one obtains the MME for the sys- 
tem's reduced density operator p{t) 



^ = -i[H,,p]+C{u;)p + C{-u)p, (2) 



where the dissipation terms are defined as follows (with 



In the above equations the following definitions are con- 
sidered: Aj{ij) = J2uj^^p I'/'a) {(t>a\ (flj + a])\(j)p) {(t)fs\ ful- 
filling the properties Aj{—uj) = A'j^uj), where U)a,i3 = 
— Oq, with Ofc as an eigenvalue of Hamiltonian Hg 
and its corresponding eigenvector \(j)k), denoting the k- 
th dressed-state. We should point out that the eigen- 
frequencics of Hamiltonian Hg are chosen in order to 
satisfy the following inequality < f^s < il4 < 
r^s < < ^1- Further in Eq. ([2]) one may use 
the so-called Kubo-Martin-Schwinger (KMS) condition 
p^ . which gives a relation for the damping constants 
7j(— w) = exp (— (Ii/Tj) 7j((I;), where Tj are the reservoir 
temperatures in the corresponding unit. The KMS con- 
dition ensures that the system tends to a thermal equi- 
librium for t — > cxD. 

In order to solve Eq. ([2]) one may use a kind of for- 
mal solution, because in the most general case there is no 
an analytic solution for the eigenvalue equation based on 
Hamiltonian Hg. Once having the operators Aj(wq,^), it 
is easy to write Eq. ([2]) for the density operator p{t) de- 
composed in the eigenstates basis, {4>m\ p{t) \4>n) = Pmm 
and we get 
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Pmn = —i^n,mPmn + ^ {'^SmeSenPkk — SmkPkn 



k=l 



— SknPmk) + ^ --^—{'2SmkSknP66 — SmdPdn — ^QnPmf^i'i) 

k=l 

Here Smn is the Kronecker 5\ the physical meaning of 
the damping coefficients 7fc_>6 and 76-).fe refers to the 
rates of the transitions between the eigenfrequencies 
rife downward and upward, respectively, defined as fol- 
lows 7fc_^6 = Z]jWi,2,3} c?7j [("('^6,fc))Tj + 1] (similar to 
Eq.(13) in Ref. |17| redefined for finite temperature) and 
76-i./£ results from the KMS condition, where Cj (with the 
integer i varying from 1 to 25) are the elements of the ma- 
trix for the transformation from the states {|1) , |6)} to 
the states {|0i) , |(/)6)} (see Eq. (14) and Appendix A 
in [13). Here (n(w„,^))T, = {e'Sip-^'^)^ _ 1)"^ corre- 
sponds to the average number of the thermal photons. 
The damping coefficients play the central role in our 
model because their dependence on the temperature of 
the reservoirs implies a complex exchange mechanism be- 
tween the elements of the system and the baths. There- 
fore, in the presence of the temperature we solve numer- 
ically the coupled system of the first-order differential 
equations and compute the evolution of entanglement 
considering the atom-field system in the initial unexcited 
state LggOOO). 

In the next section we present the calculations of the 
quantum correlations depending on the system charac- 
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FIG. 2: Evolution of the concurrence for g = u — 5^ and different atom-cavity detunings: (a) A = , (b) A = 10~'*aja and (c) 
A = O.luja- The baths have the same temperature with the average number of thermal photons given by (n(dj6,5))T- The axis 
of the dimensionless time, 'yt, is in a logarithmic scale. 



teristics, such as atom-cavity detuning, coupling con- 
stants and thermal reservoirs. In order to compute the 
atomic entanglement, we need to perform a measure- 
ment of the cavities- fiber field with a state |000) = 
|0)ci ^ |0)c2 ® \^)f- The feasibility of such a measure- 
ment is also discussed. 



III. MEASURING THE QUANTUM 
CORRELATIONS 

A. Entanglement 

Once projected on the state |000) of the field subspace 
we find that the reduced atomic density matrix in the 
two-qubit basis {|.g.g) , |.ge) , |eg) , |ee)} preserves during 
the time evolution a X-form structure 



m 



fpii ON 

p22 P23 

p32 P33 

Oy 



(4) 



with the atoms initially in ground state [i.e. pnfO) = 1]. 
The entanglement measured by the concurrence Q could 
be easily computed and gives C{t) = 2 | p23/{pii + 
P22 + P33) I ■ A more particular form of the density matrix 
(4) may result in the case of interchanging the undistin- 
guished qubits in equivalent cavities (i.e. for gi = 172, 
71 = 72 and Ti = T2). 

In the following, we are mainly interested in study- 
ing the evolution of atomic entanglement, the concur- 
rence (C), as a function of the temperatures of the ther- 
mal baths. The system under consideration refers to the 
atoms with long radiative lifetimes, each coupled to its 
own cavity. These two cavities are connected by a fiber 
with the damping rates 71 = 72 = 73 = 7 = 27r MHz, re- 
spectively, which are within the current technology [l5| . 
The transition frequency of the atom is chosen to be mid- 
infrared (MIR), i.e. LUa/^TT = 4THz and hence, for exper- 
imental purposes the coupling between the distant cavi- 
ties can be realized by using the modern resources of IR 
fiber optics, e. g. hollow glass waveguides [2^, plastic 



fibers |21|, etc. We choose the range of MIR frequen- 
cies in order to limit the thermal reservoir only up to 
room temperature (300K), which corresponds to a ther- 
mal photon. The values of the coupling constants and the 
atom-cavity detuning will be varied in order to search the 
optimal result. We must mention here that to satisfy the 
RWA we should have 2g » 'jmaxi^) [3- Satisfying this 
condition we start with the case gi = g2 = g = v ^ 67, 
considering all the reservoirs at the same temperature, 
T, and study how the atomic entanglement evolves as 
a function of the atom-cavity detuning, A. The result 
is shown in Fig. [2] from which we conclude that the 
atom-cavity detuning facilitate in this case the genera- 
tion of a quasi-stationary atomic entanglement and for 
A = O.lcja the system reaches a long-lived entanglement 
state. Of course, in the asymptotic limit the concurrence 
will vanish and the atoms eventually disentangle them- 
selves due to the damping action of the reservoirs. The 
maximal value of the concurrence of ^0.2 corresponds to 
the bath's temperature about 300 K, that is about one 
thermal excitation (consistent with the single-excitation 
approximation in the model) for the given frequency w^, 
so that kBT/hwa ~ 1.5. In order to find the optimal re- 
lation between the coupling constants and damping rate 
we did the calculations for different situations as follows: 
(i) g = y = IQO7, (ii) g = 5-f and y = IOO7, (iii) g = IO7 




FIG. 3: Concurrence for A — O.ltJa, g = 67 and v = IOO7. 
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FIG. 4: Evolution of the concurrence for arbitrary baths temperatures, (a) Ti = r2 = and varying the fiber's batli temperature, 
(b) T3 = and varying equally the cavities' bath temperatures, and (c) varying differently all the temperatures. The rest of 
the parameters are the same as in Fig. [ll^c). 



and = 7. For example, we present the case (ii) in Fig. 
[3l from which we see that the concurrence gets the same 
maximal value as in the previous case Fig. [2][c), but it 
takes a longer time for the quasi stationary entanglement 
to reach its plateau. The rest of the cases give worse re- 
sults. 



Now, let us analyze a more general situation, when all 
the independent baths have different temperatures. After 
performing the computations, we found an interesting 
effect that only the thermal bath of the fiber plays an 
important role in the generation of entanglement in the 
system, while the thermal baths of the cavities generate 
very little entanglement. This situation is represented in 
Fig. m Therefore, after analyzing all the calculations 
at the given circumstances, we come to the conclusion 
that the case represented in Fig. ^ic) corresponds to the 
optimal one for the generation of entanglement. 




yt 



FIG. 5: Evolution of the quantum discord (QD), entangle- 
ment of formation (E) and classical correlations (CC) for 
one thermal excitation and the parameters chosen as in Fig. 
[21^c). The inset represents the same quantities as a function of 
the temperatures of the reservoirs calculated for a late time, 
t = Is. 



B. Quantum Discord 

Alternatively to entanglement, the quantum correla- 
tions can be also quantified by the quantum discord 
S B H-EI ■ Since in our case the two-qubit density ma- 
trix has a simplified X form one can easily compute 
the quantum and classical correlations in the system h 
using a particular case for the algorithm discussed in 
Even if some recent studies, such as , found that the 
analytic approach of Rcf. Q could not be considered as 
a general one, in our case the computation of QD may 
follow this procedure without some divergences of the 
minimization approach. In the framework of the algo- 
rithm and notations used in [§] , we have to optimize QD 
just by changing the parameters (fc, I) in the range (0, 1) 
and found easily the condition of the resultant minimum 
for {k,l)=l/2. We have also compared the calculations 



with the approach proposed in [11 



by using Eq. (6) of 
the later and obtained exactly the same result. Hence, 
we observe in Fig. [5] the time evolution of the QD sim- 
ilar to that of entanglement, but the initial growth is 
steeper in the discord, which implies the appearance of 
the quantum correlations in the system prior to the en- 
tanglement [23 |. For a better illustration of the thermal 
effect under discussion, in the inset is shown the temper- 
ature dependence of the steady values (flat time plateau) 
of the quantum and classical correlations. 



C. Experimental hint 

In the following, we discuss the tasks important for 
an experimental realization of the ideas discussed here. 
In our opinion, the most difficult is to realize a quantum 
non-demolition (QND) measurement of the photon states 
in the fiber-coupled cavities. However, nowadays there 
exist technological possibilities to realize experiments on 
QND photon counting, attaining single-quantum resolu- 
tion, performed with optical or microwave photons (23j 
(for an exhaustive review see Ref. [131 )■ the experi- 
ment discussed in Ref. [2^ the cavity mode was coupled 
to Rydberg atoms or superconducting junctions and the 
QND method is based on the detection of the dispersive 
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FIG. 6: Probability of finding simultaneously the fiber and the 
cavities in a vacuum state by engineering of the thermal baths 
as follows: (a) all the baths have the same temperature T, (b) 
varying the fiber's bath temperature T3, while Ti = T2 = 0. 
The rest of the parameters are the same as in Fig. EJc). 



phase shift produced by the field on the wave function 
of non-resonant atoms crossing the cavity. This shift can 
be measured by atomic intcrferomctry, using the Ram- 
sey separated-oscilatory-ficld method. The advantages 
of QND experiments in radiometry and in particular ap- 
plied for IR photons are suggested in [25| . 

In order to simulate a measurement on the fiber-cavity 
subsystem one may compute the field density operator 
and therefore monitor the probability of the field state. 
As we are interested in preserving the field in the vacuum 
state, i.e. Pfib-cav{t) ~ |000) (000|, one tests the proba- 



bility of this state during the temporal evolution of the 
system. The dynamics of this probability for different 
schemes of engineering of the thermal reservoirs is shown 
in Fig. ini Based on these results we conclude that the 
success to find the fiber-cavity field in a vacuum state af- 
ter the measurement strongly depends on the managing 
of the thermal reservoirs. Hence, from this point of view, 
a more efficient variant to drive the qubits to long-lived 
quantum correlations is to increase the fiber's bath tem- 
perature while the baths of the cavities are maintained 
at the lowest possible temperature. 

IV. CONCLUDING REMARKS 

In this study we show a very interesting effect that 
the long-lived quantum correlations between the atoms 
trapped in separate cavities can be generated by the dis- 
sipative coupling to the thermal baths. This is an ex- 
ample that could give us insight into the effects of the 
system-environment exchange versus the quantum corre- 
lations. From the analysis of the obtained results, mainly 
Fig. |3]and[BJ we conclude that the entanglement can be 
optimized by engineering the thermal bath of the fiber 
rather than the baths of each cavity, hence suggesting 
that the quasilocal manipulations produce little effect on 
the generation of entanglement. Furthermore, we found 
that our system evidences quantum correlations quanti- 
fied by QD prior to the appearance of the entanglement 
(Fig. [5]). Summarizing, the model discussed here can be 
implemented as a QND measurement on the cavity-fiber 
fields with a high success probability (Fig. This is 
an example of a system where quantum correlations are 
only driven by thermal excitations and can be of interest 
as an alternative method for protection and generation 
of quantum correlations. 
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